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Abstract 

We show that a system described by an equation of state which contains a high 
number of degrees of freedom (resonances) can create a considerable amount of super- 
fluid (condensed) pions through the decay of short-lived resonances, if baryon number 
and entropy are large and the dense matter decouples from chemical equilibrium ear- 
lier than from thermal equilibrium. The system cools down faster in the presence of a 
condensate, an effect that may partially compensate the enhancement of the lifetime 
expected in the case of quark-gluon-plasma formation. 
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The investigation of hot hadronic matter through heavy ion colhsions has shown that an 
understanding of medium effects and final state interactions is of great importance for 
the interpretation of the particle spectra in order to distinguish new physical phenomena 
like the formation of quark-gluon plasma (QGP) from background phenomena caused by 
"conventional" physical effects. We refer here, e.g., to the strangeness enhancement, the 
J 1^1 suppression and the soft-pion puzzle, which can in principle be explained either by 
the formation of a QGP or from purely hadronic origin (for a recent review, see ref. Q). 

An interesting explanation for the soft pion puzzle was proposed in The authors argued 
that the baryonic and mesonic resonances may decouple from chemical equilibrium but 
still remain in thermal equilibrium. Since the resonance production rates fall drastically 
after their decoupling from chemical equilibrium, the remaining short-lived resonances 
would then rapidly decay into pions. As a consequence, the pions may acquire a non-zero 
chemical potential /Xt^ which leads to a softer p_|_-spectrum. 

A different possible explanation for the absence of local chemical equilibrium of pions in 
the final (decoupling) stage, which should not be confused with the scenario discussed in 
ref. [^, was proposed in refs.[^, ^ where it is assumed that the pions are initially created 
out of local chemical (and even thermal!) equilibrium and it turns out that the system 
never reaches complete local chemical equilibrium 

In the present letter, we investigate the scenario of ref. [^, where it is assumed that the 
dense matter reaches local thermodynamic equilibrium in an early stage of the expansion. 
The authors of [Q| considered only the case that the chemical potential stays below the 
pion mass. Moreover, they did not attempt to determine the value of the pion chemical 
potential in terms of the temperature and baryon number density at the point of decou- 
pling from local chemical equilibrium. Below, we address the problem if and under what 
conditions the decay of short-lived resonances can lead to the formation of a pionic Bose 
condensate, and what happens to the condensate in the subsequent expansion until the 
system decouples from local thermal equilibrium and the final state particles are emitted. 
To do so, we explicitly consider the transition from a resonance gas to a gas of stable and 
long-lived particles for a hadronic equation of state which contains the known resonances 
up to masses of 2 GeV, imposing constraints from the conservation of energy-momentum. 
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baryon number and strangeness. In particular, we shall show that for a system rich in 
entropy and baryon number fij^ may reach the pion mass and, consequently, the pions may 
form a Bose condensate in the central region in relativistic heavy ion collisions. 

For a hydrodynamically expanding system, the formation of such a condensate implies 
the presence of a superfluid component. The idea that the hadronic matter might be 
superfluid was put forward a long time ago in a somewhat different context Q. This 
could explain in a natural way the presence of a coherent component in multiparticle 
production. The implications for pion interferometry measurements will be discussed 
in a separate paper |^]. We note that although the creation of a remarkable amount 
of condensed pions is interesting by itself, these effects could also be used to answer 
the question of whether pions and resonances are in chemical equilibrium before they 
decouple from thermal equilibrium. Because the final stage in the hadronic phase can 
have a significant influence on the particle spectrum, the question of whether the freeze- 
out occurs in chemical and thermal equilibrium, as assumed in ^ or only in thermal 
equilibrium Q is of great importance. 

In order to investigate the conditions necessary to form a condensate, let us first consider 
a system in local thermodynamic (i.e., thermal and chemical) equilibrium. The medium is 
then completely described in terms of its equation of state (EOS) which can be expressed 
in the form 

efi = f{Tji,f,'}i) (1) 

where Sfi is the energy-density of the fluid, Tfi the temperature, and the quantities 
are the chemical potentials related to conserved charges (such as baryon number B 
and strangeness S). The index // refers to quantities which describe the fluid. 

During a heavy ion collision we expect under certain conditions the formation of a system 
of hot hadronic matter or a QGP described by eq. (||) which is in chemical and thermal 
equilibrium. The system then undergoes a phase of hydrodynamic expansion until it has 
become so dilute that the interactions are no longer strong enough to maintain local equi- 
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librium. This final stage, when the collective behaviour terminates and particles decouple 
from the dense matter, is referred to as freeze-out. 

The scenario of ref.|Q] differs from most of the other hydrodynamic descriptions of heavy 
ion collisions in so far as in the dense matter is assumed to go out of chemical equilibrium 
but remain in thermal equilibrium for some time. Or, to put it differently, there will be 
a "chemical freeze-out" independent of and prior to the "thermal freeze-out". As was 
mentioned above, the transition out of chemical equilibrium leads to the decay of short- 
lived resonances and the appearance of a non-zero chemical potential of the pions. 

Below, we shall assume that the chemical freeze-out, just like the thermal freeze-out, occurs 
"instantaneously" (i.e., on a three-dimensional hypersurface characterized by a condition 
such as T{x,t) = const.). This seems a reasonable assumption since the relevant time 
scales are the decay times of the short-lived resonances which are on the order of ~ 1 
fm/c. To illustrate the rapidness of such a transition, we consider the simplified case of a 
system of pions and rho-mesons in local thermal equilibrium that undergoes a longitudinal 
scaling ||9| expansion. In addition to the decay p — > vrvr we also take into account the inverse 
process tttt p + X. The equations which describe the time evolution of the system are 

dr ~ T 

2^-2 {av)ni (3) 



— - = — - ^ - p 

dr T 



dMp ^ _!V _ !V + ^^^^^2 

where e = £p + e-,,- and P = Pp + P-j^ are the energy density and pressure, rii {i = p, vr) 
are the number densities and r = — is the longitudinal propertime coordinate. 
Tp is the lifetime of the p, and {av) is the thermal average of the cross section for the 
process tttt p + X. We have solved eqs.(|2|- ^ numerically under the assumption that 
at T = To the vr's and p's are in chemical equilbrium at a temperature Tq, and that the 
system remains in local thermal equilibrium during the expansion (note that the eqs.(|^,^) 
imply that (n^r + 2np)T = const.). Fig. 1 shows the time evolution of the pion chemical 
potential p-,^ for tq = 1 fm/c and Tq = 200 MeV. It turns out that for reasonable values 
[av) < 100 mb, contributions of the inverse process are negligible. It can be seen that 
Ptt approaches its final value very fast, after less than 1 fm/c. For larger values of tq, the 



transition is even faster. The figure also sliows that the decay of rho- mesons alone is not 
sufficient to drive the pions into the condensate (for that, one needs to take into account 
the decay contributions of all the other short-lived resonances as well). 

We now return to the description of the expanding hot and dense matter at the point of 
chemical freeze-out. It will be assumed that the fluid freezes out into "stable" vr's , X's, 
nucleons, A's and the long living resonances w's and r/'s. The term "stable" here means 
that the decay of the particle or resonance occurs only after the complete (thermal and 
chemical) freeze-out of the system. In particular, the contribution of w's and ??'s to the 
pion chemical potential is zero. Nevertheless, their later decay leads to a non-thermal 
component of the pion spectrum which is determined by the decay kinematics [^] . 

The system of chemically frozen out particles can be described by a chemical decoupling 
temperature T'^^'f', chemical potentials for baryons and strange particles, fiB and /U5, 
and a pion chemical potential, fij^, which describes the overpopulation of pions due to the 
decoupling of resonances. The chemical potentials are determined by the requirement that 
the energy density eji, the baryonic density bji and the strangeness density sj; = of the 
fluid are equal to those of the system after chemical freeze-out. We obtain the following 
system of equations: 



.therm ( n-ich. 



-.therm/,, ,, rpch.f.\ 



+ E 



(5) 



i=K,N,K,uj,ri,... 




(6) 
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(7) 
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where 





(9) 
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are the thermal parts of the energy density and the number density of particle species i, 
and where we have used the notation 



^l^ = Y.Q^^'' (10) 

k 

for the chemical potentials. Note that on the r.h.s. of eq. (|lO|) the charges include the 
pion number A^^ which is a conserved quantity after chemical freeze-out. 

In eq. (^) we have introduced a term that describes the condensed component of pions 
which will appear if the chemical equilibrium for pions reaches the pion mass. In this case 
an overpopulation of pion states occurs and the remaining pions will be forced into the 
ground state, thereby forming a Bose condensate. The fluid then consists of a thermal 
component and a superfluid component. The latter moves with the rapidity of the fluid, 
whereas the thermal pions have on top of the fluid rapidity a thermal distribution. The 
resultant single inclusive distribution of pions emitted from a hydrodynamically expanding 
source 



1 dn ^ r pi.da'' 

P± dydpi_ (27r)3 J„ ^ p^^^^A^ ^ _ ^ 

+ J nl°''u^da''6{y - tanh (uii/uo^)6^{pi_ - m^u±) (11) 

where y and p± are the rapidity and transverse momentum of the emitted pion, is the 
four-velocity field and a is the freeze-out hypersurface. The first term on the r.h.s. of (pT]) 
is the thermal contribution given by the relativistic invariant Cooper- Frye formula ]lO[ , 
and the second term is the condensate contribution. 

The coupled system of nonlinear of equations (^)-(|7|) have to be solved in order to de- 
termine from the fluid variables eji, Ufi and sj; = the quantities ^s, and (for 
/Utt < w-tt), or and n^°" (for = m.,^) which describe the system after chemical 

freeze-out. 

Let us now investigate under what circumstances the system will develope a superfluid 
component. To this end, we apply a hadronic EOS which contains all the important 
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resonances and a treatment of compression effects, thus retaining the essential features of 



nuclear matter near the ground state 111]]. In addition to the work presented in |11| we 
take into account the effects of finite baryonic and mesonic masses, use Fermi statistics 
for the baryons and include mesons up to masses of 2 GeV. 

Fig. 2 shows the thermal and the condensate component of pions at chemical freeze-out, 
and the chemical potentials ^t^, fit, and fis, as functions of the temperature T'^^'^' that 
characterizes the decoupling from chemical equilibrium, for three different values of the 
baryon number density b. It can be seen how the baryonic and strange chemical potentials, 
/ife and /is, decrease with increasing temperature T'^^ f-. Clearly, the strange chemical 
potential induced by the finite baryon number never execeeds 150 MeV. As the amount of 
heavier resonances increases with temperature, the pion chemical potential grows with 
rpch.f.^ It becomes equal to the pion mass in the temperature region 170-200 MeV; the value 
of T'^^'f' where ^It^ = depends on the baryonic density. This behaviour is also refiected 
in the dependence of the thermal and condensate densities of pions on the temperature 
and baryon number density at chemical freeze-out. It can be seen that at sufficiently 
high baryonic densities and temperatures the condensate contribution to the chemically 
decoupled pions becomes significantly large. It is noteworthy that the condensation effect 
is enhanced if the density of the strange particles has not yet reached its equilibrium value. 
In this case the remaining energy is distributed among the non strange resonances, and 
this leads to an enhanced production of particles in the condensate. 

Finally, we need to discuss the question whether or not the condensate survives the time 
period between the chemical and the thermal freeze-out. This is of particular importance 
if the dense matter decouples from chemical equilibrium at temperatures T'^^'f' ~ 170 — 200 
MeV which are considerably higher than those usually associated with the thermal freeze- 
out, T*^-^' ~ m^r. After decoupling from chemical equilibrium, the system of stable and 
long-lived particles continues to expand until it has cooled down to the temperature T^^-f- . 

For a one-dimensional scaling expansion, the evolution of the system is determined by the 
equations 
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ST = const. (12) 

b T = const. (13) 

Se T = const. (14) 

r = const. (15) 

Ua T = const. (16) 

which describe the conservation of entropy Se, baryon number B, strangeness S and pion 
number {se, b, s and n-j^ being the corresponding densities, respectively). The index 
a labels those of the stable particle species (K, N, A,uj,rj, ...) which also have decoupled 
from chemical equilibrium. For the gas of stable particles and long-lived resonances, the 
thermodynamic quantities are given by the expressions on the r.h.s. of eqs. (|^ - ([lO|), 
with Se = J2i{^i + Pi ~ fi-i^'d/'P- Energy density and number density of the pions consist 
of a thermal and a condensate component, = e^'*^'"™ -|- rriT^n'^^ and n-,^ = n^'*^'"™ + n^°", 
respectively. Note that the condensate does not contribute to the entropy. Consequently, 
the cooling of the system - i.e., the function T(r) - does not depend on the superfluid 
component. Fig. 3 shows the cooling curves T(r) of a pion gas, for different values of 
the pion chemical potential. For comparison we have also included the results for an ideal 
gas of massless pions. It can be seen that the cooling rate increases with increasing pion 
chemical potential and becomes maximal in the presence of a condensate. This implies 
that Bose condensates may reduce the lifetime of the fireball. 

In particular, we are interested in the time dependence of the fraction of pions in the 
condensate, fcon = n'j^^ / {n^^^^"^ + n™*^). Eqs. (12,|l5) imply that for the fraction of 
thermal pions, f therm = 1 - /con, 



(nf--/se) 

(nferm/Se)\eh.f. 



rth.f. _ pch.f. \ " ' th.f. 

•I therm J therm f„therm/„ M \ ') 



where the labels th.f. and ch.f. refer to thermal and chemical freeze-out as before. The 
condensate survives until thermal decoupling if the r.h.s. of eq.([T7|) remains below 1. For 
the two limiting cases of a non-relativistic and of a ultrarelativistic gas, fconi'T) =const., 
i.e., the fraction of particles in the condensate does not change during the expansion. For 



the pion gas, it turns out that in the temperature range between 200 MeV and 150 MeV, 
ftherm Varies by less than 10%. 

To summarize, we have shown that the excited hadronic matter created in ultrarelativistic 
heavy ion colhsions could form a pion condensate if (i) the system is rich in baryon number 
and entropy density, (ii) the hadrons decouple from local chemical equilibrium earlier than 
from thermal equilibrium, as suggested in Q, and (iii) the loss of chemical equilibrium 
occurs at temperatures of ~ 180 MeV or higher. The effect depends on the chemical 
decoupling temperature, on the baryonic density and, rather sensitively, on the resonance 
contribution to the EOS.0 The effect is increased if the chemical equilibrium for strange 
particles is not complete, i.e. if the amount of strangeness which is initially zero and starts 
to increase during the collision and expansion process has not yet reached its equilibrium 
value. 

It is not the purpose of the present study to prove that such a decoupling process occurs. 
To do so would require a detailed knowledge of the density and temperature dependence 
of the inelastic and elastic cross sections of hadrons in the dense matter. Rather, we would 
like to point out here that, if this decoupling happens in a certain temperature region (like 
it was quoted in Q), then a pion condensate should appear. 

If a pion condensate is created in a heavy ion collision, it may be distinguishable from the 
case of a resonance gas in thermal and chemical equilibrium essentially by the following 
effects. 



It should lead to a coherent component which then would appear in the two-particle 
correlation functions of identical pions through a reduction of the intercept of the 
correlation function and the appearance of a two exponent behaviour of the correla- 



tion function [13|. 



The condensate component is moving with the fluid velocity and has no additional 
thermal component which smears out the distribution. This might lead to character- 



^Indeed, we have checked that for the EOS discussed in ref. [|l^, which was obtained by joining a 
parametrization of lattice QCD data to a Hagedorn type resonance gas EOS, about 40% of the pions 
would turn out to be in the condensed phase at normal freeze-out temperatures of ^ 140 MeV. 
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istic bumps and shoulders in the rapidity and transverse momentum distributions of 
pions0. One also can expect that the coherent component disappears if one considers 
particles of velocities that exceed the maximum fluid velocities. 

• The lifetime of the fireball may be reduced. This could at least partially compensate 
the enhancement of the lifetime expected in the case that quark-gluon-plasma 
is formed in the intial stage. 

Instructive and helpful discussions with R.M. Weiner, F. Navarra and H. Leutwyler are 
gratefully acknowledged. This work was supported in part by the Deutsche Forschungsge- 
meinschaft (DFG), the Federal Minister for Research and Technology (BMFT) under the 
contract no. 06MR731, the Gesellschaft fiir Schwerionenforschung (GSI) and the Los 
Alamos National Laboratory. 



■^In ref. it was suggested that at /i-^ — m-^ the thermal component may lead to a dip at px = (cf. 
also refs. given in 
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FIGURE CAPTIONS 



Figure 1: The build-up of a pion chemical potential in an expanding gas of pions and 
rho-mesons out of chemical equilibrium. 

Figure 2: Left column: Dependence of the chemical potentials for baryons (dashed), 
strangeness (dotted) and pions (solid) on the chemical freeze-out temperature T'^^-f-, for 
three different values of the baryon number density b. Right column: The thermal (solid) 
and condensate (dashed) components of the pion number densities at chemical freeze-out. 

Figure 3: The temperature evolution of a pion gas for different chemical potentials and 
condensate contributions. For comparison, the curve for an ideal gas of massless pions has 
been included as well. 
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